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\^ ■ Abstract 

£f) \ In this paper, we determine the density functions of nonsymmetrised doubly non- 

central matrix variate beta type I and II distributions. The nonsymetrised density 
functions of doubly noncentral and noncentral bimatrix variate generalised beta 
\Q • type I and II distributions are also obtained. 

o . 

q ; 1 Introduction 

When we consider generalising the distribution of a random variable to the multivariate 
case, two options are normally addressed, those of extending it to either the vectorial or the 
;_h ' matrix cases, e.g. normal, t or bessel distributions, among many others. However, some 

of these generalisations have traditionally been made directly to the matrix case, where 
such a matrix is symmetric - this is the case of the chi-square and beta distributions, for 
which the corresponding multivariate distributions are the Wishart and matrix variate beta 
distributions, respectively. Nevertheless, these latter generalisations are inappropriate in 
some cases, because sometimes we might be interested in a vectorial version and not in a 
matrix version. For example, we are int erested in a random vecto r in which each marginal 
is a random variable beta (type I or II). lLibbv and Novick ( 1982 ) proposed a multivariate 



(vector) beta dist ribution. Some app l icatio ns to utility modelling and Bay esian analysis are 



also prese nted in Libbv and Novick I (1982) and Chen and Novick (1984), respectively. In 
particular lOlkin and Liul ( 2003| ) proposed the following bivariate version. Observe that the 



following definition eliminates the hypothesis that the variables have a chi-squared distribu- 
tion, assuming, instead, a gamma distribution. 
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Let A, B and C be distributed as independent gamma random variables with parameters 
a = a,b,c, respectively and 6 = 1 in the three cases (see eq. (|4]) in Section [2]), and define 

Ul = ATc> U2 = bTc- (1) 

Clearly, U\ and U2 each have a beta type I distribution, U\ ~ BI\(a,c) and J7 2 ~ BIi(b,c), 
over < u\,U2 < 1. However, they are correlated, and then ([/i,£/ 2 )' has a bivariate 
generalised beta type I distribution over < u\, u 2 < 1. 

A similar result is obtained in the case of beta type II. Now, let us define 

A B 
F 1 --, F 2 - -. 

Once again it is easy to see that F\ and F2 each have a beta type II distribution, F\ ~ 
BII\(a,c) and F 2 ~ BIIi(b,c), over /i,/ 2 > 0. As in the beta type I case, they are 
correlated therefore and (Fi, F^f has a bivariate generalised beta type II distribution over 
A,/ 2 >0. 

These ideas can be extended to the matrix variate case. Thus, let us assume a partitioned 

matrix U = (Ui:U 2 )' € Sft 27nxm , then under the matrix variate versions of the transforma- 
tions |IJ), we are interested in finding the joint density of Ui and U 2 , where it is easy to see 
that the marginal densities of Ui and U 2 are matrix variate beta type I distributions. In 
the central case, the mat rix variate joint densities of Ui and U 2 an d of Fi and F 2 and some 
properties are studied in Diaz-Garcia and Gutierrez-Jaimezl ( 2008 ). These distributions are 



termed central bimatrix variate generalised beta type I and II distributions, respectively. 
They play a potentially important ro le in the context of shape th eory, specifically in affine 
or configuration densities, such as the lGoodall and Mardial (jl993h conjecture. Suppose that 



we have two samples of images of size n, each one of which is obtained at two times. Also, 
assume that we are interested in evaluating whether a learning process is present or whether 
the process has a memory. In this context, if we obtain as the configuration density a central, 
noncentral or double noncentral bimatrix variate generalised beta type I and II distribution, 
it might be possible to study these problems (learning or memory problems) and to compare 
the parameters of Fi (Ui) and F 2 (U 2 ) considering the latter as bimatrix variate. 

In this paper, we study bimatrix variate generalised beta type I and II distributions 
under different cases of noncentrality. Some definitions regarding the symmetrised func- 
tion are given in Subsection 12.11 and Subsection 12.21 presents known and new results about 
central, noncentral and doubly noncentral matrix variate beta type I and II distributions; 
also we include the definition of the central bimatrix variate generalised beta type I and II 
distributions. Nonsymetrised doubly noncentral density functions of the bimatrix variate 
generalised beta type I and II distributions are studied, and diverse noncentral cases of the 
bimatrix variate generalised beta type I and II distributions are obtained as particular cases 
of nonsymetrised doubly noncentral density functions, see Sections [3] and [4j respectively. 



2 Preliminary results 

2.1 Symmetrised density function 

In multivariate analysis there exist a large class of important hypothesis testing problems all 
of which may be tested by a set of criteria that depend functionally on the eigenvalues of a 
matrix var i ate. W ith the propose to investigate the non-null distributions of these criteria, 
Greenacrel (1973) introduce the notion of a symmetrised distribution of a matrix variate, 



a notion which facilitates many proofs in such derivations. 
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Given a density function / X (X), X e 5R mxm , X > 0, iGreenacrd || 1973ft proposes the 
following definition 



/.(X)=/ / x (HXH')(rfH), HeO(m) 



where O(m) 



'O(m) 

{H e j R mxm |H / H = HH' 



I m } and (dH) denotes the normalised Haar 



measure on 0(m), see iMuirheadl (|1982l pp. 60 and 260). This function / S (X) is termed 
symmetrised density function of X. 

Our proposal is to apply this idea from Greenacrel (1973) in an inverse way, i.e. well- 
known the explicit expression of the symmetrised density function of X 



/.(X) = / /(HXH')(dH). 

JO{m) 



(2) 



We wish to identify the density function /(X). The density function obtained by applying 
the idea underlying $2$ is termed the nonsymmetrised density function. Finally, note that the 
joint density function of the eigenvalues of X can be found from / S (X) or /(X), indifferently. 



2.2 Matrix variate beta distributions 



In general, matrix variate beta type I and II distributions are defined in ter ms of two matri- 
ces, sa y, A and B, whi c h are independent a nd have W isha rt distributions , s eelOlkin and Rubin 
1964ft . [Khatril Jl970h. IMuirheadl (Il982h. [Farrelll (Il985h. ICaded (| 1996ft . iGupta and Nagarl 
2000), iDiaz-Garcia and Gutierrez-Jaimea (|2007l 120061 120081 ). among many others. The 



present paper generalises these results, assuming that A and B have matrix variate gamma 
distributions. 

The m x m matrix A is said to have a noncentral matrix variate gamma distribution 
with parameters a € 3? in which is an m x m positive definite matrix and H is an m x m 
matrix , this fact being den oted as A ~ Q m (a : 8, SI) , if its density function is (see Muirheadl 
(|l982l . pp. 57 and 61) and lGupta and Nagarl (|200Cf )) 



g m (A; a, 0, O) = g m (A; a, &) Q F 1 (a, SI® 1 A), A > 0, 



(3) 



where o-Fi( 4 ) is a hypergeometric function with a matrix argument (see Muirhead ( 19821 p. 
258)) and Q m (A;a, 0) = Q m (A;a, 0, 0) denotes the density function of a central matrix 
variate gamma distribution given by 



| A |a-(m+l)/2 

g m (A;a,0) = ' 1 etr(-0- 1 A), A > 0, 

I m [aJ|0| Q 



(4) 



and denoted as A ~ Q m {a, 0) = G m {a, 0, 0). Where etr(-) 
the multivariate gamma function and is defined as 



exp(tr(-)) and T m [a] denotes 



T m [a] = f etr(-V)|Vr-( m+1 )/ 2 (rfV), 

JV>0 



Re(o) > (m- l)/2. 

In addition to the c lassi fication of beta distributions as beta type I and type II (see 

iGuot a and Nagarl (|2000l ) andlSrivastava and Khatrl t wo de fi nitions have been pro- 

posed for e ach of these, se e Olkin and Rubinl (| 1964ft . Srivastava ( 1968ft , Diaz-Garcia and Gutierrez- Jaimez 
( 2001 ) and Jamesl (1964). Let us focus initially on the beta type I distribution; if A and B 
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have a matrix variate gamma distribution, i.e. A ~ Q m (a,I m ) and B ~ Q m (b, I m ) indepen- 
dently, then the beta matrix U can be defined as 



U = 



(A + B)- 1 /2 A (( A + B)- 1 / 2 )', Definition 1 or, 
AV2( A + B)" 1 (A 1 /2)' j Definition 2, 



(5) 



where C 1 / 2 (C 1 ^ 2 ) / = C is a reasonab le n qnsingular f actori zation of C, see lGupta and Nagar 
(2000), ISrivastava and Khatr i (1979) and iMuir head! l|l982l ). It is readily apparent that under 
definition 1 and 2 the density function is 



BI m (V;a,b) 



1 



Pm [a, b] 



| U |a-(m + l)/2| Im _ U| b -( m+1 )/ 2 (dU), < U < I„ 



(6) 



writing this fact as U ~ BI m (a,b), with Re(a) > (m — l)/2 and Re(6) > (m — l)/2; where 
/3 m [a, 6] denotes the multivariate beta function defined by 



Pm[b, a] 



| S |a-(m+l)/2| Im _ g|6-(m+l)/2( dS J 



0<S<I„ 



| R |°-(m+l)/2| Iri 
R>0 

r m [a]r m [b] 

T m [a + b] 



■R|-( Q+fc )(dR) 



A similar situation arises with the beta type II distribution, and thus we have the fol- 
lowing two definitions: 



B- 1 / 2 A(B- 1 / 2 )', Definition 1, 
AVSb-^A 1 / 2 )', Definition 2, 



(7) 



with the distribution being denoted as F ~ BII m {a,b). In this case, under definitions 1 and 
2, the density function of F is 



BII m (F;a,b) 



1 



| F |o-(m+l)/2| I 



-(a+b) 



F > 0. 



(8) 



/3 m [a,b] 

Diaz-Garcia and Gutierrez- Jairne3 (2007. 2006) showed that in doubly noncentral and 



noncentral matrix variate beta type I and II distributions, the corresponding density func- 
tions are invariant under definitions 1 and 2. Therefore, henceforth we shall make no dis- 
tinction between definitions 1 and 2. 

When these ideas are extended to the doubly noncentral case, i.e. when A ~ Q m {a, I m , Oi) 
and B ~ G m (b, I m , O2), strictly speaking, we have not found the densities of the matrix vari- 
ate beta t ype I and II distributions. Rather, for the case of the beta type II distribution, 
(Chikusd . Il980h found the distribution of V = B- 1 / 2 A(B- 1 / 2 )' where A = H'AH and 
B = H'BH, H 6 0(m), with 0(m) = {H e Sft m xm |HH / = H ; H = It is straight- 

forwar d to show that the procedure proposed by I Chikusd ( 1980f) and Chikuse and Davisl 
( 1986) is equival ent to fi nding the symmetrised density de fi ned b y iGreenacrd (|l973f ). see 
alsolRouxj (|l975l ). From iDiaz-Garcfa and Gutierrez- Jaimezl (2006) and using the notation 
for the operator sum as in David ( 19801 ) we have the following: 

1. The symmetrised density function of doubly noncentral matrix variate beta type I is 

BJ m (U;a,6)etr(-(fii+f! 2 )) (9) 



{a + b) q 



OO 

k £ {a) K {b) x k\ n 



C<p (Im ) 



< U < L, 
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2. and the symmetrised density function of doubly noncentral matrix variate beta type 
II is 

BJJ m (F;a,6)etr(-(ni+n 2 )) (10) 



E 

K..A; 



(a + b)* C£' A (ni, n 2 )C£' A ((I m + F)-*F, (I m + F)- 1 ) 



(a) K (6) A fc! I! 



C0(I m ) 



where F > 0, Re(a) > (to — l)/2, Re(fe) > (m — l)/2, (a) T is the generalised hypergeo- 
metric coefficient or product of Pochhammer sy mbols and C^ ' A (-,-) d enotes th e inva riant 
polynomials with the mat rix arguments defined in Davis (1980), see also Chikuse (1980) and 
Chikuse and Davis! (jl986t) . 



As particular cases of doubly noncentral distributions it is possible to obtain two different 
definitions of noncentral distributions, giv en an other classification , in wh ich the beta matrix 
is defined as follows, sec Grecnacrc (|l973h and lGupta and Nagari (|2000h : 

W = A 1/2 (A + B)- 1 (A 1/2 )', denoting as BI(A) m (a,b,n 2 ), («i = 0) 
U = A 1/2 (A + B) -1 (A 1/2 )'> denoting as BI(B) m (a, b, n x ), (fi 2 = 0). 

Similarly, in the case of beta type II we have 

V = B- 1/2 A(B~ 1/2 )', denoting as BII(A) m (a, b, Sl 2 ), («i = 0) 
F = B- 1 / 2 A(B- 1 / 2 )', denoting as BII(B) n (a, b, Hi), (fl a = 0) 

Both distributions, types A and B, play a fun damental role in various areas of stati stics, for 
example in the W and U criteria proposed by iDiaz-Garcia and Caro-Lopera ( 20081 ) . 

The symmetrised and nonsymmetrised density functions of W, U, V and F can be ob- 
tained as particular cases of © and (fTU)) . All these densities are found in lDiaz-Garcia and Gutierrez- Jaimez 
(|2007f ). 

Now, using the approach described in IDiaz-Garcia and Gutierrez- Jaimed (j2007t ) we can 
find the (nonsymmetrised) density functions of do ubly noncentral ma trix variate beta type 
I and II distributions. Observe that for ([9]) and bv lDiaz-Gareial(|2006l ). 



C2' A (niHUH', n 2 (I m - HUH'))(dH) 



O(m) 



where 6^' A is defined in lDavid (|l979h and lChikusd (|l980h . 

Proceeding in analogous form for (|10[) , we have the following. 



Theorem 2.1. for Re(a) > (to - l)/2 and Re(6) > (to - l)/2, 

1. the nonsymmetrised density junction of the doubly noncentral matrix variate beta type 
I is 

BI m (\J;a,b)etT(-(n 1 + n 2 )) (11) 



E 



(a + 6), 



c^ x (n x u, n 2 (i m - u)) , o < u < i r , 



(a)«(6) A fc! J! * 



which is denoted as U~ BI m (a, b, fli, fl 2 ). 

I?, and i/ie nonsymmetrised density function of the doubly noncentral matrix variate beta 
type II is 



£// m (F;a,&)etrHft 1 + n 2 )) 



(12) 
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c«' A (n 1 (i ro + F)- 1 F,n 2 (i r , 



F > 0, 



which is denoted as U ~ BII m (a, b, S~2 2 ). 



Doubly noncentral, noncentral and central matrix variate beta type I and II distribu- 
tions play a very important role in diverse problems for proving hypotheses in the context of 
multivariate analysis, including canonical correlation analysis, the g eneral linear hypo t hesis 
i n MANOVA and multip le matrix variate corr elation analysis, see iMuirhead (1982), Rao 
( 1973 ). Srivastava ( 19681 ) and iKshirsagar ( 196lh . Similarly, doubly noncentral and noncen- 
tral beta distributions are to b e found in the contex t of ec onometrics and shape theory, see 
Chikuse and David (Il986l) andlGoodall and Mardial (Il993l) . respectively. 



Now from Diaz-Garcia and Gutierrez- JaimezT i 20081 ): let A, B and C be independent, 
where A - G m (a,I m ), B - Q m (b,I m ) and C ~ G m (c,I m ) with Re(a) > (m - l)/2, Re(6) > 
(to — l)/2 and Re(c) > (to — l)/2 and let us define 

Ui = (A + C)" 1/2 A(A + C)" 1/2 and U 2 = (B + C)" 1/2 B(B + C)" 1/2 (13) 

Of course, Ui ~ BI m (a, c) and U2 ~ BI m (b, c). However, they are correlated and therefore 

the distribution of U = (Ui:U2)' G gfj 2mxm [ s termed a central bimatrix variate generalised 
beta type I distribution, denoted as U ~ BQBl% m -^ m {a^ b, c). Moreover, its density function 



is 



| Ui |a-(m+l)/2| U2 |b-(m+l)/2|j m _ Ui |6+c-(m+l)/2 Ij^ _ jj 2 |o+c-(m+l)/2 



U!U 2 



(3^[a,b,c]\I, n 

and is denoted as ,8S£>/ 2mxm (U; a, b, c), where < Ui < I„ 
(to - l)/2, Re(6) > (m - l)/2 and Re(c) > (m - l)/2 and 



(14) 



< U 2 < I m with Re(a) > 



/3 r *„[a,6,c] = 



r m [a]r m [b]r m [c 

T m [a + b + c] 



Similarly, let 



Clearly, Fi 



Fi = C-^AC- 1 / 2 and F 2 = C-^BC- 1 / 2 (15) 
BII m (a,c) and F 2 ~ BII m (b,c). But they are correlated and then the 



distribution of F = (Fi:F 2 )' G 5R 2mxm is termed a central bimatrix variate generalised beta 
type II distribution, which is denoted as F ^ BQBIl2mxm(a, b, c). Its density function is 



BGBIh m xm{^a,b,c) 



| Fl |o-(m+l)/2|p 2 |6-(m+l)/2 

/3£ l [a,&,c]|I tn + Fi+F 2 |°+H-° 



(16) 



where Fi > 0, F 2 > 0, with Re(a) > (to - l)/2, Re(6) > (to - l)/2 and Re(c) > (to - l)/2. 

Other properties of bimatrix variate generalise d beta type I and II distributions are 
studied in Dia z-Garcia and Gutierrez - Jaime d (|2008l ). 

The use of matrix and bimatrix variate beta-type distributions has not been developed 
as expected and hoped, due particularly to the fact that such distributions depend on hy- 
pergeometric functions with a matrix argument or on zonal polynomials, which until very 
recently were quite complicated to evaluate. Recently, descriptions have been made of al- 
gorithms that are very efficient at calculating both zonal polynomials and hypergeometric 
functions with a matrix argument; th ese can be us ed more widely and more efficie ntly in 
noncentral distributions in general, see iKoev (2009) and IKoev and Edelman I (|2006l ). 
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3 Doubly noncentral bimatrix variate generalised beta 
type I distribution 

In this section we derive the doubly noncentral bimatrix variate generalised beta type I 
distribution. 

Theorem 3.1. Let A, B andC be independent random matrices, suchthat A ~ Q m (a, l m ,Sli), 
B ~ G m (b, I m , r2 2 ) and C ~ G m (c, I m , tt 3 ) with Re(a) > (m — l)/2, Re(6) > (m — l)/2 and 
Re(c) > (m — l)/2 and Zei ms define 

Ui = (A + C)- 1/2 A(A + C)- 1/2 and U 2 = (B + C)" 1/2 B(B + C)" 1/2 (17) 



T/ien t/ie symmetrised density function oft] = (Ui:U 2 )' £ 5ft 2 



fl0BJ aro xm(U;o 1 6 > c)etr{-(ni + n 3 + n 3 )} ]T -J^±^ - (18) 

c;- T - A (ni, n 2 , » 3 )g;- T ' A (M 1 , m 2 , m) 

CV(I m ) 

and then the nonsymetrised density function o/U = (Ui:U 2 )' G 3fJ 2mxm is 

iBaB/ 2n , XTO (U;o ) 6,c)etr{-(n 1 + n 2 +n3)} V ,° ... , C * * , n ( 19 ) 

x Cj'^COxMi.naMa.naM), 
which is denoted as U ~ £>C?£>i 2mX m(a, b, c, fii, 1~2 2 , $13); where < Ui < I TO and < U 2 < 
I m and 

Mi = (Im-UaJ^-UiUa)- 1 ^, 
M 2 = (I rn -U 1 )(I m -U 1 Ua)- 1 U 2) 
M = (^-Ui^-UiUaJ-^^-Ua). 

Proof. The joint density function of A, B and C is 

ctr{-(n 1 + » 2 + »3)} _ (m+1)/2 b _ (m+1)/2 _ (m+1)/2 + B + C)} 

r m [ajr m [6jr m [cj 

x 0^(0; fiiA) i?i(&; O2B) Fi(c; n 3 C), 
Now, consider the transformations (fT3f and C = C, then 



(20) 



(dA)(dB)(dC) = |C| m + 1 |I m -U 1 r( m + 1 )|I m -U 2 |-( m + 1 )(dU 1 )(dU 2 )(dC). 
The joint density function of Ui, U 2 and C is 

etr{-(fli + Vt 2 + » 3 )} ixj^a-Cm+i^i^-Cm+i)^ 

r m [a]r m [6]r m [ c ] |i m - u^p^+W^; _ u^p^+ip 

x |C| a + h + c -( m+1 )/ 2 exp{-(I m - Ux)- 1 ^ - U!U 2 )(I m - Ua)"^}, 
x 0F1 (a; fliC 1 / 2 ^ _ Uxj-^iC 1 / 2 ) oi^i (&; OaC 1 ^^ - U 2 )- 1 U 2 C 1 /2 

x Fi (c; n 3 C) 

Now, note that 
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n _tt 'i-i-T -/ I 1 ™ _ u i) 



i) 



With similarly expressions for (I m — U2) 1 — I m - 

ii) From the argument of etr(-) in pi] ), denotes M^ 1 = (I m -Ui) _1 Ui + (I m -U 2 ) _1 U2+I I1 
and note that 



1VT 1 



(I m - Ui)" 1 ^ - UiU 2 )(I m - Us)- 1 , 
(I m - Ua)" 1 ^ - U!U 2 )(I m - Ui)" 1 . 



iii) Now assuming that (I m — Ui) 1 UiM is an argument of a symmetric function (/(AB) 
/(BA)) by i) and ii) we have 



Ui^UiM 



' (I m -V2)(J m -ViV2)- 1 Vi 
(Im-UiUaJ-^Im-UajUx, 
UiCIn,-U 2 )(I m -UiU 2 )- 1 , 

UiCI n ,-UiU 2 )- 1 CI ff ,-U 2 ). 



With similar expressions for (I m — U 2 ) 1 U2M. 
The marginal joint density function of Ui and U 2 is 

etr{-(»i + » 2 + »3)}|Ui| a ~ (m+1)/2 |U2| fc ~ (m+1)/2 

r m [a]r m [&]r m [ C ]|I m - _ U 2 |&+(m+l)/2 

x / |C| a+frfc - (m+1 )/ a eBcp{-M- 1 C}oi ; i (ojniC^f^-UiJ^UiC 1 / 2 ), 
ioo ^ ' 

x F! (6; fiaC 1 / 2 ^ - U 2 )- 1 U 2 C 1 / 2 ) ^i (c; ft 3 C) (dC), 

Denoting the density joint function of Ui and U 2 as f v U2 (Ui, U 2 ), considering the corre- 
sponding symmetrised function 

/.(U 1> U 3 )= / / „ j (HUiH' I HU 2 H')(dB) 

JO(m) 

and the transformation C = HCH' with (dC) = (d(HCH')). Then expanding the hyper- 
geometric functions o-Fi(') in terms of zonal polynomials we obtain that / s (Ui,Ua) is 

etr{-(Oi + n 2 + 3 )}|U 1 | Q -( m + 1 )/ 2 |U 2 | fc -( m + 1 )/ 2 



r m [a]r m [6]r m [ c ]|i m - Ui|°+(-+ 1 )/ 2 |i m - u 2 | b +(-+ 1 )/ 2 

1 



00 00 00 



k=0 t=0 1=0 k t A 



{a) K (b) T (c) x k\t\l\ 



I |c| Q + b + c -( m+1 )/ 2 exp{-M- 1 C}, 
Jc>o 



I C K ffirHC 1 / 2 ^ - UO^UiCVH') , 

Jo(m) V ' 



x C T (flaHC 1 / 2 ^ - U^UaC^Hj C x (n 3 HCH') (dH)] (dC), 
By integrating with respect to H, using Chikus e and David (|l986l . equation (2.2)) and the 
notation for the operator sum as in iDavid (|l980h~ we have 

etr{-(ni + n 2 + n 3 )}|Ui| a -( m+1 )/ 2 |U2| Mm+1)/2 



r m [o]r m [6]r m [ c ]|i m - u x |«+(-+i)/2|i m - u 2 |"+(™+d/2 
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CT X ((I™ - U^-iUxM- 1 , (I m - U 2 )- 1 U 2 M- 1 , M- 1 ) 

x q^ v (dC) - 



Finality, by integrating with respect to C, see IChikusd (| 19801 equation (3.21)) and iii), we 
obtain the joint symmetrised density function of Ui and U2. 

The joint nonsy mmetrised density function of Ui and U2 is obtained by applying the idea 
of lGreenacrei(|l973h in an inv erse way. With this p roposal observe that BQBl2mxm (HUH'; a, b, c) — 
B(/<B7 2TOXm (U; a, b, c) and bv lDiaz-Garcial (|2006l ) 

f C/ r ' A (rJiHMiH', rJiHM 2 H', niHMH')(dH) 

JO(rn) 

_ C£ T ' x (Sl u n 2 , » 3 )C7' A (M 1; M 2 , M) 
0% T ' C^(I m ) 

from where the desired result is obtained. □ 
In addition, note that in Theorem 13. 1[ 



Ui ~ BI m (a, c, ni,0 3 ) and U 2 ~ BI m (b, c, fi 2 , fi 3 )- 

Next, assuming that one and/or two of the matrices A, B or C have a central matrix 
variate gamma distribution in Theorem 13 - 1 1 let us study all the possible nonsymetrised 
noncentral densities. 



Corollary 3.1. Let us assume the hypothesis of Theorem \3.1\ Then the joint nonsymetrised 
density function of Ui and U 2 is: 

1. if Ui = n 2 = 

BgBl2mx m (U;a,b,c) 

1-F1 (a + b + c;c; S2 3 M) 

etr{S2 3 } 

and Ui - BI(A) m (a, c, n 3 ) and U 2 - BI(A) m (b, c, to 3 ). 

2. Iffl 3 = 

BgBI 2m xm(U;a,b,c) ^ (a + 6 + c) <% T 

etr{n 1 + n 2 } \ (aUb) T k\t\ C * ^M^M,) 

also we have that XJ\ ~ BI(B) m (a, c, J7i) and U 2 ~ BI(B) m (b, c, J7 2 ). 

3. 7/ fi 2 = o 3 = 

B£B7 2TOXm (U;a,6,c) A/T x 

iFi(a + + c; c; SZiMi) 

etr{S2i} 

and with Ui ~ BI(B) m (a, c, Oi) and U 2 ~ BI m (b, c). 

4. // Hi = n 3 = 

S£B7 2TOXm (U;a,&,c) OM > 

1F1 a + 6 + c; c; 2 M 2 ) 

etr{S2 2 } 

and Ui ~ BI m (a, c) and U 2 ~ BI(B) m (b, c, fi 2 ). 
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5. if n 2 = o 

In addition, note that Ui ~ BI m (a, c, ill , O3) and U2 ~ BI(A) m (b, c, O3). 

6. I/fii = 

etr{n 2 + n 3 } 2. ( 6 ) T ( c)A t!Z! C * ^M 2 ,n 3 M) 

and where Ui ~ BI(A) m (a, c, Jls) and U2 ~ BI m (b, c, O2, n 3 ). 

Proof. The joint density functions of Ui and U 2 in al l items are a consequence of the 
basi c propert i es of invariant polynomials, see iDavisl ( 19791 equations (2.1) and (2.3)), see 
also Chikusei ( 1980l equations (3.3) and (3.6)). The second claim in each of the items is a 
consequence of construction (flTj) . □ 



4 Doubly noncentral bimatrix variate generalised beta 
type II distribution 

In this section we derive the doubly noncentral and noncentral bimatrix variate generalised 
beta type II distributions. 

Theorem 4.1. Let A., B andC be independent random matrices, such that A ~ Q m (a,l m ,Sli), 
B ~ G m (b, I m , r?2) and C ~ Q m {c, I m , ^3) with Re(a) > (m — l)/2, Re(6) > (m — l)/2 and 
Re(c) > (m — l)/2 and fei us define 

Fj = C-^AC" 1 / 2 and F 2 = C-'^BC" 1 / 2 (21) 



T/ien i/ie symmetrised density function o/F = (Fi:F 2 )' S sft2mxm j s 

BffB/I amxm (P;a,6,c)etr{-(n 1 + n a + n,)} £ n 0®) 

k,t,A;0 

cE' r,A (ni, n 2 , n 3 )C2 ,T ' A (Ni, n 2) n) 



C<^> (Im ) 

and £/ien the nonsymetrised density function o/F = (FxiF2)' G sft 2mxm j s 

BgB/I 2fflXm (F;a > 6 > c)etr{-(n 1 + n 2 + n 8 )} V ^ ... . „ (23) 

(a) K (o) T (c)AK! *! /! 

x ^ T ' A (f2 1 N 1 ,02N 2! 3 N), 
which is denoted as F ~ BGBIl2mxm(a, b, c, Ox, $12, O3); where Fx > and F2 > and 

Nx = (I m +Fx)- 1 (I m + F2)(I m + F 1 +F 2 )- 1 Fi, 
N 2 = (I m +Fi+F 2 )- 1 (I m + F 1 )(I m +F 2 )- 1 F2, 
N = (I m +F 1 +F 2 )- 1 . 
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Proof. The proof is similar to that given for Theorem 14.11 Or alternatively, by noting 
that if U ~ BI m (a,b), then (I m — U) -1 — I m ~ BII m (a, 6), see lSrivastava and Khatril (|1979l ) 
and Diaz-Garcia and Gutierrez- Jaimez ( 20071 ). Then, by construction (|17p we have that if 
U ~ ^/ 2mxm (a,6,c,ni,n 2 ,0 3 ) then F = (F 1 |F 2 )' ~ B0B// 2TOXro (a, b, c, fl u fi 2 , n 3 ), 
where 

Fi \ / (I m - U0- 1 - 1„ 
F 2 J I (^-U,)- 1 -I„ 



F 



Then the inverse transformation is given by 



Ui \ / I m — (I m + Fi) 1 



U 2 / 1 J ™ - i 1 ™ + F 



2; 



1 



(24) 



and the Jacobian is given by 

(dUi)(dU 2 ) = |I m +F 1 |-( ro+1 )|I m + F 2 |-( m+1 )(dF 1 )(dF 2 ). 
Then, (f2"2")l follows, observing that under transformation (flM)) 
i) 

ggg/2mxm((Im + (Im + Fi)- 1 ^ - (I m + Fa)' 1 )';^ c) 

|I ro -Fi|('»+i)|I ro +Fa|("*+ 1 ) 

= BgBII 2m xmiV;a,b,c), 

Mi = (I m + F 1 )- 1 (I™+F 2 )(I m + F 1 +F 2 )- 1 F 1 = N ls 
M 2 = (I m + F 1 +F 2 )- 1 (I m +F 1 )(I m + F 2 )" 1 F 2 = N 2 , 
M = (I m + F 1 +F 2 )- 1 = N. 



ii) and 



To obtain (f2"5|) , observe that 

BQBIhm x m (HFH' ; a, b, c) = BSS// 2mXm (F; a, b, c) 
and bv lDiaz-Gareial (|2006f ) 

/ C"' T ' A (fiiHNiH', fJ!HN 2 H', n x HNH')(dH) 

JO(m) 

_ c^isi^ n 2 , » 3 )g;' T ' A (N 1 , N 2 , N) 

from where the nonsymmetrised joint density function of Fi and F 2 is obtained. □ 

Finally, let us find the different possibilities of the nonsymetrised noncentral density 
functions of bimatrix variate generalised beta type II distributions. 

Corollary 4.1. Under the hypothesis of Theorem \4-l\ the joint nonsymetrised density 
function of Fi and F 2 is: 



1. if fix = n 2 = 

BGBII 2mxm {¥;a,b,c) 



x F x (a + &+c; c; £l 3 N) 



etr{ft 3 } 

andFi ~BII(A) m (a,c,n 3 ) andF 2 ~ BII(A) m (b,c,n 3 ). 
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2. if n 3 = o 



BgBII 2mxm (¥-a,b,c) ™ (a + b + c)^ 6 K * 

2^ "Twirrnr^ (S2iJNi,s2 2 jn 2 ) 



ctr{Oi + n 2 } 



(a) K (b) T k\ t\ * 



In addition, we have that Fi ~ BII(B) rn (a, c, fii) and F 2 ~ BII(B) m (b, c, ^2). 
5. 7/ 2 = 3 = 

BgBII 2mxm {¥;a,b,c) 

^(a + 6 + c; c; O lNl ) 

and wi£a Fi ~ BII(B) m (a, c, Jli) and F 2 ~ BII m (b, c). 

4. if rii = n 3 = 



S£S77 2mXm (F;a,&, c) 



iFi(a + 6 + c; c; 2 N 2 ) 



ctr{0 2 } 

and Fi - BII m (a, c) and F 2 - BII(B) m (b, c,fl 2 ). 

5. If 2 = 

ggg77 2mxm (F;a,b ; c) ^ (a + M^ , 

7n addition note that Fi ~ BII m (a, c, fii, O3) and F 2 ~ BII(A) m (b, c, ^3). 

6. 7/fii = 



BGBII 2mxm (¥;a : b,c) ^(a + b + c)^ T / 



E 



CT' A (n 2 N 2 ,n 3 N) 



etr{ft 2 + tt 3 } ^ (6)t(c)a*!« 
and where Fi ~ BII(A) m (a, c, $^3) and F 2 ~ BII m (b, c, i~2 2 , $13). 
Proof. The proof is obtained in a similar way to the proof of Corollary 13.1 



□ 



5 Conclusions 

The problem to finding the density function of a d oubly n oncen tr al beta t ype I I dis- 
tribution has been studied by di f ferent authors, see Davisl ( 19791 ). Chikuse (|l980l) and 
Diaz-Garcia and Gutierrez- Jaimezl (2006), among others. All these authors, in fact, found 
the symmetrised density function of a doubly noncentral beta type II distribution. The 
nonsymmetrised density function remained an open problem. The Theorem 12.11 solves this 
problem for doubly noncentral beta type I and II distri butions by app lying in an inverse 
way the definition of s ymmetrised function proposed by IGreenacrd ll 19731) . This approach 
was used previously by iDiaz-Garcia and Gutierrez- Jaimezl ( 2007 ) in the noncentral cases. 

In a similar way, we have found the symmetrised doubly noncentral generalised beta 
type I and II distribu tions and, by ap plying in an inverse way the definition of symmetrised 
function proposed by IGreenacrd ( 19731) . we have found the noncymmetrised doubly noncen- 
tral generalised beta type I and II distributions. Finally, as corollaries, we studied all the 
different noncentral generalised beta type I and II distributions. 
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